Optimally robust shortcuts to population inversion 
in two-level quantum systems 



(N 



O 

C 
3 



A. Ruschhaupt 1 ' 2 , Xi Chen 34 , D. Alonso 5 and J. G. Muga 3 ' 4 

1 Institut fur Thcoretische Physik, Leibniz Universitat Hannover, Appelstr. 2, 30167 
Hannover, Germany 

2 Department of Physics, University College Cork, Ireland 

3 Departamento de Qufmica-Ffsica, UPV-EHU, Apdo 644, 48080 Bilbao, Spain 
QO 1 4 Department of Physics, Shanghai University, 200444 Shanghai, China 

5 Departamento de Fi'sica Fundamental y Experimental, Electronica y Sistemas and 
C~| ' IUdEA, Universidad de La Laguna, 38203 La Laguna, Spain 

Abstract. We examine the stability versus different types of perturbations of 
recently proposed shortcuts-to-adiabaticity to speed up the population inversion of 
a two-level quantum system. We find optimally robust processes using invariant based 
engineering of the Hamiltonian. Amplitude noise and systematic errors require different 
optimal protocols. 
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\Q ! 1. Introduction 

o 

Manipulating the state of a quantum system with time-dependent interacting fields is a 
fundamental operation in atomic and molecular physics, with applications such as laser- 
controlled chemical reactions, metrology, interferometry, nuclear magnetic resonance 
(NMR), or quantum information processing [TJ |2j [31 H] . For two-level systems there are 
several approaches proposed to attain a complete population transfer, for example, n 
pulses, composite pulses, adiabatic passage and its variants. In general, the n pulses 
may be fast but highly sensitive to variations in the pulse area, and to inhomogeneities 
in the sample pp. Used first in nuclear magnetic resonance [5], composite pulses provide 
an alternative to the single 7r-pulse, with some successful applications [6J [7], but still 
need an accurate control of pulse phase and intensity. A robust option is in principle 
adiabatic (slow) passage, which is however prone to decoherence because of the effect 
of noise over the long times required. A compromise is to use speeded-up "shortcuts to 
adiabaticity" , which may be broadly defined as the processes that lead to the same final 
populations than the adiabatic approach in a shorter time. 

Several methods to find shortcuts to adiabaticity have been put forward [8j El QUI 
[TTl [T2| [T3| [HI [15j [16j [TTJ, [18] for two- and three-level atomic systems. The transitionless 
or counter-diabatic control protocols, proposed by Demirplak, Rice [8] and Berry [9] 
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start from a reference time-dependent Hamiltonian H and provide an extra interaction 
that cancels the diabatic couplings. This results in an exact following of the adiabatic 
dynamics of the reference Hamiltonian, in principle in an arbitrarily short time. They 
have been applied, for example, to speed up the RAP for an Allen-Eberly scheme 
[TO] . Modified by a unitary transformation [16], the transitionless quantum driving 
has been experimentally implemented for a two-level system realized by Bose-Einstein 
condensates in optical lattices [T5] . 

Another shortcut technique is to inverse engineer the Hamiltonian using Lewis- 



Riesenfeld invariants [19], as in [2Ql [2H [221 E31 [2H [25J [261 EZl [28l [29l [30l [3H [32]. 



The invariant-based method has been applied to accelerate the adiabatic processes 
for trap expansion or compressions [20l [ZD, (221 [231 1211 (2S1 [26] and atomic transport 
[2?1 EH [29] . It has also been combined with optimal control theory [23j [29] , and proposed 
for other applications [301 EH E21 [16]. Counterdiabatic and invariant-based engineering 
can in fact be shown to be potentially equivalent methods by properly adjusting the 
reference Hamiltonian [13J. In standard applications though, Hq is set according to 
some predetermined, standard protocol (for example Landau-Zener, Allen-Eberly, or 
finite-time schemes), and the formulation and results of the two methods are generally 
quite different, so they may be considered in practice separate approaches. 

A key element to choose among the fast protocols is their stability or robustness 
versus different perturbations. We will compare the results with ordinary (flat) n 
pulses and explore the stability of the transitionless approach with respect to parameter 
variations for a finite-time sinusoidal protocol for Hq. 

The main aim of this paper is to find optimal protocols with respect to amplitude 
noise of the interaction and with respect to systematic errors. The optimality will 
be determined by minimizing properly defined sensitivities. It turns out that the 
perturbations due to noise and systematic errors require different optimal protocols, 
and we shall use invariant-based inverse engineering to find them. 

The rest of the paper is organized as follows. In the following section we shall review 
the transitionless-based shortcuts protocol and the invariant-based one. In section [31 the 
general formalism to model amplitude-noise error and systematic error will be presented. 
The special case of solely amplitude- noise error will be examined in section [4] where the 
noise sensitivity of the different protocols will be studied and the most stable protocol 
will be derived. In section O the special case of solely systematic error will be studied 
and the most stable protocol will be derived. The general case of amplitude-noise as 
well as systematic noise will be for the different protocols will be numerical studied in 
section [6l 

2. Shortcuts to adiabatic passage for a two-level quantum system 

We assume a two-level system with a Hamiltonian of the form 



H (t) 



h I -A(t) 
2 \ Sl R (t) + iSlj(t) 



n R {t) -in T (t) 

A(t) 



) 



(i) 
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For example, in quantum optics such a Hamiltonian describes the semiclassical coupling 
of two atomic levels with a laser in a laser-adapted interaction picture. In that setting 
Q(t) = Qn(t) +ifli(t) would be the complex Rabi frequency (where Qr and Qi and the 
real and imaginary parts) and A would be the time-dependent detuning between laser 
and transition frequencies. We find it convenient to keep the language of the atom- 
laser interaction hereafter noting that in other two- level systems Q(t) and A(t) will 
correspond to different physical quantities and that instead of "atom" one may refer, 
for example, to a spin-1/2, or to a Bose-Einstein condensate on an accelerated optical 
lattice [15] . 

Initially at time t — 0, the atom is in the ground state. Often the goal is to achieve 
a perfect population inversion such that at a time t = T the atom should be in the 
excited state. The time T should be as small as possible but also the scheme or protocol 
to achieve this population inversion should be as stable as possible concerning errors. 
In the following subsections we will review different schemes to achieve a population 
inversion before we discuss different types of possible error sources in the next section. 

2.1. tt pulse 

A simple scheme to achieve population inversion is a ir pulse. In this case the laser is on 
resonance, i.e. the detuning is zero A(£) = for all t. If the Rabi frequency is chosen 
like fl(t) = \£l(t) \ e ia , with a time-independent a, and such that 



the population is inverted at time T. A simple example is the "flat" 7r pulse with 



2.2. Adiabatic schemes and transitionless shortcuts to adiabaticity 

The population inversion may also be achieved by an adiabatic scheme. Let the 
instantaneous eigenstates of the Hamiltonian Hq be \n(t)), with n — 0, 1. The adiabatic 
theorem tells us that if we start in an eigenstate at t — 0, i.e. |^(0)) = \n(0)) 
and if we vary the Hamiltonian infinitesimally slowly, then the system will stay in 
the corresponding instantaneous eigenstate for all times, up to a phase factor, i.e. 
\ip(t)) ~ e lKn ^\n{t)) . If the eigenstate corresponds initially to the ground state and 
at t = T to the excited state (up to a phase) then we would achieve a perfect population 
inversion as T -> oo. 

Demirplak and Rice [8] and independently Berry [9] proposed a modification of the 
Hamiltonian such that the state would exactly follow the instantaneous eigenstate of the 
Hamiltonian Hq for an arbitrary duration T. If the desired time evolution operator is 
U — y~l e lKn ®\n(t))(n(0)\, the corresponding Hamiltonian leading to this time evolution 



is H 0a (t) = ih(d t U)U+. We may write H 0a = H + H a , where H a = ihJ2 n \d t n{t))(n{t)\. 
is the "counter-diabatic" (CD) term that guarantees that the system will follow the 




(2) 



7rt gia 
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instantaneous eigenstates of H without transitions even for a small T. This method is 
thus termed counter-diabatic approach or transitionless-tracking algorithm. 

For the two-level system with Hamiltonian Ho and Qi = this additional 
Hamiltonian takes the the form 

H(t)- h ( ° - iQa{t) \ (i\ 
with Q a = [QrA — QrA]/Q 2 . The total Hamiltonian is therefore [10] 

^ ( V*" ) • w 

which we will call transitionless shortcut protocol in the following. 
2.3. Inverse engineering of invariant-based shortcuts 

Shortcuts to adiabaticity can be also found making explicit use of Lewis-Riesenfeld 
invariants. For the general Hamiltonian H in (pQ), a dynamical invariant of the 
corresponding Schrodinger equation (this is a Hermitian Operator I(t) fulfilling J^I + 
1[Hq, I] = 0, so that its expectation values remain constant) is given by 

h f cos (0(t)) sin (9(t)) W \ 
l j 2^ ^ sin(9(t))e-W -cos(6(t)) )' { > 

where \i is an arbitrary constant with units of frequency to keep lit) with dimensions 
of energy, and the functions Q(t) and a(t) satisfy the differential equations 

6 = Qj cos a — Qr sin a, , . 

a = — A(i) — cot (cos a + sin a Vli) . 
The eigenvectors of the invariant are 

| cos (6/2) e~ ia l 2 \ ., ( sin (0/2) e~ ia l 2 , 

'^>=( sm(e/2)eW 2 )' = ( -a»(0/2)^ • ' (?) 



with the eigenvalues A general solution of the Schrodinger equation can 

be written as a linear combination = c + e tK+ ^ \<j> + (t)) + c_e* K_( ^|0_(t)), where c± 

are complex, constant coefficients, and k± are the Lewis-Riesenfeld phases [19] 

1/ - 1 u '+/■ 



d 

ih- - h 

^4 - Ho 



1 

In particular we may construct the solution 

-nW/2 _ ( cos(6/2)e— / 2 ^ 
sin(e/2)e 

and the orthogonal solution (for all times (i/)±\ip) = 0) 



W)) = \Mt))e"^ = [ TAo',r.,„/ 2 I «- ,l/2 (8) 
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where 

7 = —2n + = 2k_. 
Finally, we get 

7 = — - (cos a Qr + sin a fij) . (10) 

sin 6 

Equivalently we may design a solution of the Schrodinger equation \tp(t)) with the 
parameterization of a pure state given in (jSJ). (Note that \ip(t))(ip(t)\ is a dynamical 
invariant.) By putting this ansatz into the Schrodinger equation, we get immediately 
(J6]) and (TTOl) . A solution which is orthogonal to (151) , i.e. (V>±|V ; ) = for all times, is then 
directly given by Q). 

The next step to find invariant-based shortcuts is to inverse engineer the 
Hamiltonian. For achieving a population inversion, the boundary values should be 
0(0) = and 6(T) = tt, so 

\m) = ( e "7 /2 ) e_l7(0)/2 ' ^ = ( ) e ~ n(T)/2 - 

Assume that 0(t) and are given. Then we get for the Hamiltonian corresponding 
to this solution by inverting ([6]). This leads to 

Qr = cos a sin© 7 — sin a 0, (11) 
fli = sin a sin 9 7 + cos a 6, (12) 
A = — cos0 7 — a. (13) 

By implementing these functions exactly the population would be inverted in the 
unperturbed, error-free case. Note that invariant-based shortcuts and transitionless 
shortcuts may be formally related, see [T3] . 



3. General formalism for systematic and amplitude-noise errors 

We shall now consider systematic errors as well as noise-related errors. Let the ideal, 
unperturbed Hamiltonian be H . For systematic errors, the actual, experimentally 
implemented Hamiltonian is if 01 — Ho + 0Hi, but the evolution of the pure quantum 
state is still described by the Schrodinger equation with the perturbed Hamiltonian 
Hqi. Sometimes systematic errors cannot be avoided, for example if different atoms 
at different positions are subjected to slightly different fields, due to, for example, the 
Gaussian shape of the laser inducing different Rabi frequencies. It is thus desirable to 
have protocols which are very stable with respect to perturbed Hamiltonian functions. 

The second type of error is a stochastic one, i.e. the Hamiltonian H m is perturbed 
by some stochastic part XH 2 describing amplitude noise. A stochastic Schrodinger 
equation (in the Stratonovich sense) is then 
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dWt 



,, is heuristically the time-derivative of the Brownian motion W t . We 



where 

have (£(£)) = and (£(£)£(£')) = 5(t — t') because the noise should have zero mean 
and the noise at different times should be uncorrelated. If we average over different 
realizations and define p(t) = (p^) then p(t) satisfies 



!p = -i[ff«,d-£[/W*„A 

More details on the derivation can be found in the appendix. 

We may consider the two effects together with the master equation 



(14) 



d 

dt p 



% -[H Q + pH u p\-^[H 2 ,[H 2lP \\, 



(15) 



where /3 is the amplitude of the systematic noise described by the Hamiltonian Hi and 
A is the strength of the amplitude noise. 

In this paper, we assume that the errors affect the frequencies fl R and Qi but not 
the detuning A, which, for an atom-laser realization of the two-level system is more 
easily controlled. For the systematic error we restrict ourselves to an error Hamiltonian 
of the form 



h 



H (t) 



A=0 



(16) 



o n R (t) - zfij(t) 
n R (t) + isii(t) o 

For the noise error we restrict ourselves to independent amplitude-noise in Q R as well 
as in Qj with the same intensity A 2 , i.e. the final master equation is 



7 >2 

Jt P = - l j[H, + /?#!, p}-— 2 ([H 2R , [H 2R , p}} + [H 2I , [H 2I , p}}) 
where 



(17) 



H 2R {t) 



o n R (t) 
n R (t) o 



, H 2I (t) 



-^(t) 
iQ T (t) 



A motivation for this modeling is that two different lasers may be used to implement 
the two parts of the Rabi frequency. 

It is now convenient to represent the density matrix p(t) by the Bloch vector 



f(t) 



' Pl2 + P21 

i{pn ~ Pti) 

\ P U - P22 



\ 



J 



such that p = |(1 + r ■ a) where a = (o~i, o~ 2 , 0-3) are the Pauli matrices. The Bloch 

(18) 



equation corresponding to the master equation (117D is 

L + pLi - X 2 L 



d _ 
dt 



where 



Lr 



( A{t) 
-A(t) 

^ -Qj(t) n R (t) 



-n R (t) 
/ 



/ 



Li 













y -Q T {t) n R (t) 



-n R (t) 
/ 
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and 





V 



o 

n R (ty 



o 
o 



\ 



o n R (t) 2 + nj(t) 2 J 



Note that the probability to be in the excited state at time t is P 2 (t) = |(1 — r 3 (t)). In 
the following section we will first study the amplitude- noise errors only, then in section 
the systematic errors and finally both together. 



4. Amplitude-noise errors 



We assume that there is an amplitude-noise type or error affecting the Rabi frequencies 
and no systematic errors (/3 = 0). Let us define the noise sensitivity as 



Qn 



1 d 2 P 2 



2 d\ 2 



dP 2 



A=0 



where P 2 is the probability to be in the excited state at final time T, i.e. P 2 ~ 1 — g^A 2 . 
A smaller value of means less sensitivity with respect to amplitude-noise errors, i.e. 
the scheme is more stable concerning this type of noise. In general an analytic solution 
of the master equation ffTTl) or the Bloch equation (IT51) cannot be found. To calculate q^ 
we do a perturbation approximation of the solution keeping only terms up to A 2 (with 
/3 = 0). In this manner we get 

r 3 = (0,0, l)f 

^r 0i3 (T)-A 2 / dt'(0,0,l)U (t,t')L 2 (t')f (t') 
Jo 

^r 0i3 (T) + A 2 / dt'(0,0,-l)U (t,t')L 2 (t')r (t'). 
Jo 

where Uq is the unperturbed time evolution operator for the Bloch vector. If the noiseless 
scheme works perfectly, i.e. r ^(T) = — 1, then 

\2 l-T 
P 2 = 1 - - 

2 

where r means the transpose operation and the noise sensitivity becomes 



dt'r (t') T L 2 (t')r (t'), 



1 

~ 4 



dtf (t) T L 2 (t)f (t) 



dt 



+l\2\ 



(19) 



4-1. Example: ir pulse with real Rabi frequency 

As a first simple example of a population-inversion protocol we look at a it pulse with 
a real Rabi frequency, i.e. we set A = 0, fli = and J Q Qn(t)dt = ir. In this 
case the master equation resp. the Bloch equation for amplitude noise can be solved 
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Pi 




Figure 1. (Color online) Probability P2 to end in the excited state at time T versus 
noise parameter A. Optimal protocol (green, dashed line), flat 7r-pulse with purely 
real Rabi frequency (blue, dotted line), pure adiabatic (black, dashed-dotted line), 
transitionlcss shortcut method (red, solid line). Additional parameter for adiabatic 
protocol and transitionless protocol: f2 T = 5.57/4.3 n, SqT = (5.57/4.3) 2 7r. 



analytically. The solutions of this equation with the initial conditions pn(0) = 1, 
Pi2(0) = P2i(0) = P22(0) = resp. r(0) = (0, 0, l) r at initial time t = are 

n(t) = 

r 2 (t) = - e- x2 ft <*<* , >*7* sin na{t ' )dt ' 
r 3 (t) = e -A 2 /o T ^(*')^72 cos f£n a (f)dt?j , (20) 



which yields 



The noise sensitivity is now 

I r T 

II -2 



q N = - I Q 2 R (t')dt', (21) 



2 - 

which may be bounded as |j= < < \ max < t < T |f2 R (t)|, where the lower bound is 

2 

J dtfl R (t) < J dtfl 2 R (t). We can achieve the 



derived using Schwartz inequality, i.e. 



lower bound using a constant Q R = it /T (i.e. a flat 7r-pulse). The excitation probability 
P 2 for this flat 7r-pulse is plotted in figure [H versus the noise intensity A (blue, dotted 
line). The noise sensitivity is q^ = tc 2 /4:T 2.467/T. The other lines in figure [1] 
correspond to different protocols, see below for more details. The important thing at 
this point is to note that the stability of a protocol is very well quantified by qN, which 
is the curvature at A = 0. 

4-2. Example of a transitionless protocol 



As another example, we will now look at the stability and noise sensitivity of a 
transitionless shortcut protocol. Our reference scheme is the finite-time sinusoidal model 
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EES [31 



= H sin 



(22) 



f J , A(t) = -(jocos^ 

with £7/ = and < t < T. The excitation probability is also shown in figure [T] (black, 
dashed-dotted line). The chosen intensities are not large enough and the population 
inversion is not complete. 

Figure [T] shows the excitation probability also for the transitionless shortcut based 
on this sinusoidal model (red, solid line). The noise sensitivity of this protocol is 
q jv = 3.21/T. Figure |2] shows the noise sensitivity for the transitionless protocol based on 
the sinusoidal model ( )22|) for different values of 5$ and Qq. The minimal noise sensitivity 
in this figure is achieved for Sq = Q$ = 0.5 /T and it has the value = 2.475/T which 
is very close to the noise sensitivity of the flat n pulse in the previous subsection. 



4-3. Optimal scheme 

We can also write the unperturbed Bloch vector in the form 

/ sin 6 cos a \ 



r (t) 



V 



sin sin a 
cos O 



(23) 



This Bloch vector corresponds to the pure state in ([8]). Therefore we get from the Bloch 
equation the same equations ([6]). If the trajectory of the Bloch vector r(t) and A(t) is 
given, i.e. 9, a and A are given, then the corresponding Qr and Qi can be calculated 
by (HU) and flTJJ. Let ra(t) = tan0(A + a). Using (jl9"]h and (fTT])-(fl3lh we get for 
the noise error sensitivity 

-T 



qN 



dt 



[cos + cos a sin 0)(msina — cosa;0) z 



+ (cos 2 + sin 2 a sin 2 0)(m cos a + sin a©)" 



dtL(m, a, 9, 9), 



(24) 
(25) 
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where L is the Lagrange function for q^. We are looking for functions m(t), 0(i), ct(t) 
which minimize this functional. From the Euler-Lagrange formalism we get 

dL sin(4a) sin 2 

= — — =>- m - 



Moreover 



dm 4 cos 2 + 2 sin 2 (2a) sin 2 ' 



9L 

= — sin(4a) = =>- a = mr/A. 



da 

From this it also follows that m(t) = 0. Finally we have 

dL d dL , . 

0= ae-JtW (26) 

Let us now consider the cases n odd and n even separately. 

Case n even If n is even, then (T26]) simplifies to = 0. Taking the boundary conditions 
0(0) = 0, 0(T) = tt into account, we arrive at 

0(t) = id /T. 

It follows that 

/ 77,7r \ tt fmr\ TT 

iiR = — sin — — , Ut — cos — — , A = 0. 
V4/T V4/T 

Note that either Q,r or Qj is zero, so these schemes are flat tt pulses with a purely real 

or purely imaginary Rabi frequency. (As an example, we get for n = 6 a tt pulse with 

a flat, real Rabi frequency = ^ and Qj = 0.) For all these schemes an analytical 

solution of the master equation can be derived similar to the one in subsection 14.11 and 

the noise sensitivity of all schemes is q^ = tt 2 / (4T). 



Case n odd For n odd we get 

(3 + cos(20))0 = sin(20)(0) 2 (27) 

Then Vt R = ±Q/y/2 = ±flj and A(t) = 0. We first solve (jZZD f o r @ numerically and then 
put the solution in the expression for the noise sensitivity. The numerically calculated 
0(t) can be seen in figure [3] (solid line, left axis). The corresponding Rabi frequencies 
for n = 7 are f^(t) = Qi(t) = Q(t), where Q is shown also in figure [3] (dashed line, right 
axis). Note that for other values of odd n only the signs of Qr resp. Qt are switched. 
The noise sensitivity value is q^ = 1.82424/T < 7r 2 /(4T). Therefore, for n odd, smaller 
noise sensitivities can be achieved than for n even, so these protocols are least sensitive 
to amplitude noise. Finally, the optimal tt pulse is shown (green, dashed line) in figure [fl 
it has a noise sensitivity q^ = 1.82424/T. 

Note that an approximate solution of ( 1271) is given by 0(t) = irt/T — sin(27rt/T), 
with a noise sensitivity of q^ = 1.82538/T. 
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Figure 3. (Color online) Protocol with minimal noise sensitivity q^, @(t) (blue, solid 
line; left axis), f2(i) (red, dashed line; right axis). 



5. Systematic errors 



In this section we shall only consider systematic errors, i.e. A = 0. It is enough to work 
with pure states, instead of density matrices, satisfying the Schrodinger equation 

ihj t \i>(t)) = (H (t)+pH 1 )\i(>(t)). 
We define the systematic error sensitivity as 



1 d 2 P, 



qs 



2 <9/3 2 



0=0 



0=0 d(P 2 

where P 2 is the probability to be in the excited state at final time T. 
Using perturbation theory up to 0((3 2 ) we get 



dtC/o(T,t)tfi(t')|^o(t)> 



where \ipo{t)) is the unperturbed solution and (Jq the unperturbed time evolution 
operator. We assume that the error-free (j3 = 0) scheme works perfectly, i.e. \ipo(T)) = 
e 8/i |2) with some real [i. Then, 



p 2 = i(2|^(T))i 2 = wnMnmtmn 



1 - 



t 

h 2 



T 2 

dt(Mt)\Hi(t)\Mt)) 



because U (s,t) = \ipo(s)) (tp (t)\ + \4>±(s)){ip±(t)\, where (4 , ±it)\ip {t)) = for all times 
and \if)j_(t)) ^ s a ^ so a solution of the Schrodinger equation, see (j9]). From this we get the 
systematic-error sensitivity value 



qs = h 



T 



dt{^{t)\H y {t)\^{t)) 
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Figure 4. (Color online) Excitation probability P2 versus systematic-error parameter 
/3: protocol with zero systematic-error sensitivity (blue, dashed-dotted line), 
transitionless protocol (red, solid line), n pulse with minimal noise sensitivity (green, 
dashed line). 

5.1. Example: n pulse 

Let A = 0, Sl(t) = \£l(t)\e ia and / Q alt = tt, which correspond to a 7r pulse. Then 

we get that Hq = Hi and an analytical solution exists, 

P2 = \ ~ \ cos ((1 + 0) jf \Q(t')\dt^j = I - i cos ((1 + /3)tt) . 

It follows that qs = 7r 2 /4 independently of the time duration T. 

The excitation probability versus systematic noise is shown in figure HJ As an 
example we are looking at the ir pulse which was optimal for amplitude-noise error 
in the previous section (green, dashed line). It has the systematic-error sensitivity 

2 

qs = \ ~ 2.47 which is equal for all tt pulses. Even if the protocol is maximally robust 
concerning amplitude-noise, it is very sensitive to systematic errors. 

5.2. Example of a transitionless shortcut 

We again look at the example of a transitionless shortcut based on the sinusoidal model 
which was examined in subsection 14.21 The excitation probability versus systematic 
noise /3 is shown in figure H] (red solid line). The transitionless shortcut based on the 
sinusoidal model is more stable concerning systematic errors that any it pulse. 

Figure [5] shows the systematic-error sensitivity for the transitionless-based protocol 
for different values of S and fl . Again, the protocol takes as a reference the sinusoidal 
model (|22l) . Note that the systematic-error sensitivity qs for any it pulse corresponds 
to the upper x-y plane in this figure. This means that for all parameters shown the 
transitionless shortcut is less sensitive (i.e. more stable) concerning systematic error 
than any it pulse. 
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Figure 5. Systematic-error sensitivity qs versus fio and So for the transitionless 
protocol; the systematic-error sensitivity for a ir pulse corresponds to the upper x-y 
plane. 



5.3. Optimal scheme 



To find an optimal scheme we shall use the invariant based technique. The pure state 
\ip{t)) can be parameterized as in (|S]). The boundary values should be 0(0) = and 
O(T) = 7T. We get for the functions in the Hamiltonian leading to this solution 

Qr = cos a sin 6 7 — sin a 
Qj = sin a sin 6 7 + cos a 
A = — cos 6 7 — a 

Note that, contrary to section HI it is now more convenient to take j(t) as a given 
function instead of A(t). A solution which is orthogonal to ([H]), i.e. = for all 

times, is given by fl9]). The expression for the systematic error sensitivity is now 



Qs 



dt{* ± (t)\H 1 (t)/h\ip(t)) 



dt 



dt 



-ie * 7 7 cos 9 sin 6 + e n 6 
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dt 



; C osesine) + e- i7 e 



where we have applied partial integration in the last step taking into account the 
boundary values 0(0) = and 0(T) = ir. The expression can be further simplified 
and we get finally 

dte'^Q sin 2 . 



Qs 



independently of 



In the special case when 7(t) is constant in time, we get qs = 
0(t). With the choice a(t) constant we recover the 7r pulse. 

The minimum of qs is clearly achieved if qs = 0. In the following we will show 
that there are protocols which fulfill this condition, i.e. protocols maximally stable with 
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Figure 6. (Color online) Rabi frequency Sl/j (red, solid line) and detuning A (blue, 
dashed line) for a protocol which has zero systematic-error sensitivity. 



respect to systematic errors. We will give an example class which fulfills qs = 0. Let 

7(t) =n(2e-sin(26)). 
For this choice of 7 we get 
sin 2 (mr) 

qs = ^^- 

So for n — 1, 2, 3, ... we get protocols fulfilling q$ = 0. Note that in the limit of n — > 
(i.e. 7 — > 0), we get qs — >■ which is consistent with the previous paragraph. The 
functions in the Hamiltonian in this case are 

VLr = (An cos a sin 3 6 — sin a) O, 
flj = (4n sin a sin 3 + cos a) 0, 
A = — 4n cos sin 2 — a. 

Note that this class of protocols might not be the only ones fulfilling qs = 0. There 
is still some freedom left. For example, one could in addition require that A = or 
Qi = 0. In the following we will look at the second condition, i.e. Qi = for all t. This 
leads to a(t) = — arccot (4nsin 3 0). In addition, there is the freedom to chose 0(t) with 
0(0) = 0, 0(t/) = 7r. For n — 1 and 0(i) = 7rt/t/, the resulting Rabi frequency and the 
detuning are shown in figure |6l 

6. Systematic and amplitude-noise errors 

Finally, we will consider both types of errors together. Optimal schemes in this case 
would depend on the ratio between amplitude-noise error and systematic error in the 
experiment. We will just examine numerically the behavior of some protocols with 
respect to amplitude-noise and systematic error. Specifically we compare the minimal 
noise error protocol, the minimal systematic error protocol, and the example of a 
transitionless shortcut studied before, see figure [3 The figure shows that the different 
optimal schemes perform better than the other one depending on the dominance of one 
or the other type of error. 
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Figure 7. (Color online) Probability P2 versus noise error and systematic error 
parameter; (a) transitionless protocol (red), optimal systematic stability protocol 
(blue), optimal noise protocol (green); same result as contour plots: (b) transitionless 
protocol, (c) optimal systematic stability protocol, (d) optimal noise stability protocol. 



Summarizing, in this paper we have examined the stability of different fast protocols 
for exciting a two-level system with respect to amplitude-noise error and systematic 
errors. First we have looked at the noise error alone and we have introduced a noise 
sensitivity. We have shown that a special type of tt pulse is the optimal protocol with 
minimal noise sensitivity. Then we have looked at the systematic error alone and we 
have introduced a systematic error sensitivity. We have shown that there are protocols 
for which this sensitivity is exactly zero. Finally, we have looked at the general case 
with noise and systematic errors together. 

Future work may involve extending the present results to different types of noise 
and perturbations. The existence of a set of optimal solutions for systematic errors 
also opens the way to further optimization with respect to other variables of physical 
interest. 

Appendix A. Derivation of the Master equation for Amplitude-Noise Error 

The evolution of the quantum state with amplitude noise can only be described by a 
master equation [35]. We assume that the Hamiltonian has a deterministic part H m 
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and a stochastic part containing XH 2 . We need a mapping from a fixed time to another 
infinitesimally close, so our starting point will be 

where dt is the infinitesimal time step and dWt the corresponding noise increment in 
the Ito sense. The properties of such noise are: (dW) = 0, (dW 2 ) = dt. If we expand in 
Taylor series (lA.ljl and keep terms up to first order in dt and dW (using the Ito calculus 
rules) we arrive at the following Stochastic Schrddinger equation (SSE) 

\ 2 ' \ 

m = ~H 01 dt\i>) - 2fi2#fd#> - jH 2 dW t \ij). (A.2) 



The master equation derived from this SSE is then (JTJ 

An equivalent approach in the Stratonovich sense is to start from 

/ \ 

dW t 



dt 

V IS J 

where £(t) is heuristically the time-derivative of the Brownian motion W t . We have 
= an d = ~~ because the noise should have zero mean 

and the noise at different times should be uncorrelated. If we average over different 
realizations and define p(t) = (p%) then p(t) is fulfilling (|T4j) . To show this we define 
Pt(t) = We start from the dynamical equation for namely 

|ft = -ip«,«l-^.f«], (A.3) 
that after averaging over the noise becomes 

j t P = ~[H up}-j[H 2 ,(^)}- (A-4) 
Novikov's theorem applied to white noise takes the form 

1 / 5F 

mm) 



2 \6£(s) 
Using it we get 



= P], (A.5) 

which leads to (JHJ). 
Acknowledgments 

We acknowledge funding by Projects No. GIU07/40, No. FIS2009-12773-C02-01, No. 
FIS2010-19998, NSFC No. 61176118, and the UPV/EHU under program UFI 11/55. 



[1] Allen L and Eberly J H 1987 Optical Resonance and Two-level Atoms (New York: Dover) 

[2] Vitanov N V, Halfmann T, Shore B W and Bergmann K 2001 Annu. Rev. Phys. Chem. 52 763 

[3] Bergmann K, Theuer H and Shore B W 1998 Rev. Mod. Phys. 70 1003 



Optimally robust shortcuts to population inversion in two-level quantum systems 17 



[4] Krai P, Thanopulos I and Shapiro M 2007 Rev. Mod. Phys. 79 53 
[5] Levitt M 1986 Prog. Nucl. Magn. Reson. Spectrosc. 18 61 

[6] Collin E, Ithier G, Aassime A, Joyez P, Vion D and Esteve D 2004 Phys. Rev. Lett. 93 157005 
[7] Torosov B T, Guerin S and Vitanov N V 2011 Phys. Rev. Lett. 106 233001 

[8] Demirplak M and Rice S A 2003 J. Phys. Chem. A 107 9937; 2005 J. Phys. Chem. B 109 6838; 

2008 J. Chem. Phys. 129 154111 
[9] Berry M V 2009 J. Phys. A 42 365303 
[10] Chen X, Lizuain I, Ruschhaupt A, Guery-Odelin D and Muga J G 2010 Phys. Rev. Lett. 105 
123003 

[11] Masuda S and Nakamura K 2010 Proc. R. Soc. A 466 1135; 2011 Phys. Rev. A 84 043434 
[12] Muga J G, Chen X, Ibahez S, Lizuain I and Ruschhaupt A 2010 J. Phys. B: At. Mol. Opt. Phys. 
43 085509 

[13] Chen X, Torrontegui E and Muga J G 2011 Phys. Rev. A 83 062116 

[14] Ibahez S, Martmez-Garaot S, Chen X, Torrontegui E and Muga J G 2011 Phys. Rev. A 84 023415 
[15] Bason M G, Viteau M, Malossi N, Huillery P, Arimondo E, Ciampini D, Fazio R, Giovannetti V, 

Mannella R and Morsch O 2012 Nat. Phys. 8 147 
[16] Ibahez S, Chen X, Torrontegui E, Ruschhaupt A and Muga 2011 J G larXiv:1112.5522l 
[17] Fasihi M A, Wan Y D and Nakahara M 2011 arXiv.1110.6707. 
[18] Lacour X, Guerin S and Jauslin H R 2008 Phys. Rev. A 78 033417 
[19] Lewis H R and Riesenfeld W B 1969 J. Math. Phys. 10 1458 

[20] Muga J G, Chen X, Ruschhaupt A and Guery-Odelin D 2009 J. Phys. B: At. Mol. Opt. Phys. 42 
241001 

[21] Chen X, Ruschhaupt A, Schmidt S, del Campo A, Guery-Odelin D and Muga J G 2010 Phys. Rev. 

Lett. 104 063002 
[22] Chen X and Muga J G 2010 Phys. Rev. A 82 053403 
[23] Stenfanatos D, Ruths H, Li J S 2010 Phys. Rev. A 82 063422 

[24] Schaff J F, X.-L. Song, P. Vignolo, and G. Labeyrie, Phys. Rev. A 82, 033430 (2010). 
[25] Schaff J F, Song X L, Capuzzi P, Vignolo P and Labeyrie G 2011 EPL 93 23001 
[26] Schaff J F, Capuzzi P, Labeyrie G and Vignolo P 2011 New J. Phys. 13 113017 
[27] Torrontegui E, Ibahez S, Chen X, Ruschhaupt A, Guery-Odelin D and J. G. Muga J G 2011 Phys. 
Rev. A 83 013415 

[28] Torrontegui E, Chen X, Modugno M, Schmidt S, Ruschhaupt A and Muga J G 2012 New J. Phys. 
14 013031 

[29] Chen X, Torrontegui E, Stefanatos D, Li J S and Muga J G 2011 Phys. Rev. A 84 043415 

[30] Li Y, Wu L A and Wang Z D 2011 Phys. Rev. A 83 043804 

[31] del Campo A 2011 Phys. Rev. A 84 031606(R); 2011 EPL 96 60005 

[32] Choi S, Onofrio R and Sundaram B 2011 Phys. Rev. A 84 051601 (R) 

[33] Lu T, X. Miao X, and H. Mctcalf H 2007 Phys. Rev. A 75 063422; Lu T S 2011 Phys. Rev. A 84 
033411 

[34] Miao X, Wertz E, Cohen M G, and H. Metcalf H 2007 Phys. Rev. A 75 011402(R) 
[35] Carmichael H J 1999 Statistical methods in quantum optics 1 (Springer) 



